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The ’Brillouin zone spin filtering’ mechanism of enhanced tunneling magnetoresistance (TMR) is 
described for magnetic tunnel junctions (MTJ) and studied on an example of the MTJ with licp Co 
electrodes and hexagonal BN (h-BN) spacer. Our calculations based on local density approximation 
of density functional theory (LDA-DFT) for Co(0001)/h-BN/Co(0001) MTJ predict high TMR 
in this device due to Brillouin zone filtering mechanism. Owning to the specific complex band 
structure of the h-BN the spin-dependent tunneling conductance of the system is ultra-sensitive 
to small variations of the Fermi energy position inside the BN band gap. Doping of the BN and, 
consequentially, changing the Fermi energy position could lead to variation of the TMR by several 
orders of magnitude. We show also that taking into account correlation effects on beyond DFT 
level is required to accurately describe position of the Fermi level and thus transport propertied 
of the system. Our study suggests that new MTJ based on hep Co-Pt or Co-Pd disordered alloy 
electrodes and p-doped hexagonal BN spacer is a promising candidate for the spin-transfer torque 
magnetoresistive random-access memory (STT-MRAM). 

PACS numbers: 73.40.Rw, 85.75.-d 


Theoretical prediction of high TMR in Fe/MgO/Fe 
MTJ due to so-called ’symmetry spin filtering’ mecha¬ 
nism 0, 0 and its quick experimental verification [i, 0 
revolutionized the hard disk drive (HDD) industry during 
the last decade. On the other hand, continues progress in 
new areas of magnetic memory technology (for example, 
STT-MRAM) demands for a development of novel MTJ 
in which high value of the volume type magnetocrystalline 
anisotropy (MCA) of electrodes is as important as high 
TMR. Standard nowadays Fe-based electrodes have low 
volume MCA due to high (cubic) symmetry and thus has 
to rely on interface anisotropy which cannot fully satisfy 
strict demands of new technology. 

In present paper we suggest different mechanism of 
high TMR in magnetic tunnel junctions which we call 
the ’Brillouin zone filtering’ and study this mechanism in 
details on a specific example of low symmetry (hexago¬ 
nal) hcp-Co/h-BN/licp-Co junction. Specifically, we de¬ 
fine the ’Brillouin zone filtering’ as a mechanism of en¬ 
hanced TMR in electrode/spacer/electrode tunneling de¬ 
vice where spin filtering is provided by (1) existence of 
so-called ’hot spot’ - special part in the in-plane two- 
dimensional Brillouin zone (2D BZ) where the semicon¬ 
ductor spacer has very high probability of transmission, 
and (2) absence of states with in-plane momentums fc 
corresponding to the ’hot spot’ in one spin channel of 
the ferromagnetic electrode and presence of states with 
corresponding fc|| in another spin channel (see Fig 1 for 
illustration). Since finite area of the 2D BZ is involved, 
the Brillouin zone spin filtering mechanism potentially 
leads to exponential increase of TMR with semiconductor 
thickness, t , TMR oc exp(2A r yt), where Ay is the differ¬ 
ence of the minimal attenuation constants of the spacer 




FIG. 1: (color online). Illustration of the ’Brillouin zone spin 
filtering’ mechanism of enhanced TMR (see text). 


in areas of the 2D BZ available for majority and minority 
electrons of the electrode. This is much stronger depen¬ 
dence compared to linear increase of TMR with semicon¬ 
ductor thickness originated from symmetry spin filtering 
where filtering mechanism works only in close vicinity of 
some high-symmetry points of 2D BZ (e.g. T point in 
conventional Fe/MgO/Fe MTJs 01 )- 

While the BZ filtering is general mechanism of en¬ 
hanced spin dependent transport, in present letter we 
focus on novel physics of realization of the mechanism 
in tunneling devices. We show that the BZ filtering con¬ 
ditions can be satisfied in p-doped BN MTJ with hep 
Co electrodes. The many-electron perturbation theory 
was used to quantify excited electron states and thus find 
conditions for spin dependent transport enhanced beyond 
that of reported for CoFe/MgO system. In addition, low 
crystal symmetry realization of BZ filtering mechanism is 
of significant potential technological interest in the con¬ 
text of novel STT-MRAM technology that has a potential 
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to become an ’universal memory’ Q combining all the 
strengths and none of the weaknesses of existing memory 
types. 

In STT-MRAM devices the perpendicular magnetic 
anisotropy (PMA) of electrodes is preferable option as 
compared to electrodes with in-plane anisotropy due to 
faster switching with low current, higher thermal sta¬ 
bility and scalability @0- The requirements on PMA 
electrodes include high thermal stability at reduced di¬ 
mensions, low switching current and high TMR all at the 
same time. Various PMA materials have been studied for 
STT-MRAM electrodes, including (Co,Fe)/(Pt,Pd) mul¬ 
tilayers @-0, Ll 0 -ordered (Co,Fe)Pt allo ys loUTT! . rare- 
earth/transitional-metal (RE/TM) alloys [12,[l3j: ultra- 
thin CoFeB 0 and CoFeAl jliiHa] films, and tetragonal 
manganese alloys such as M^Ga 17]. Despite of inten¬ 
sive research none of these materials satisfy strict require¬ 
ments that would allow STT-MRAM to replace conven¬ 
tional memory today. Electrodes with large concentration 
of heavy Pt or Pd elements (like multilayers Co/Pt or 
Co/Pd) have high PMA but also exhibit relatively large 
Gilbert damping constant due to strong spin-orbit cou¬ 
pling of Pt and Pd, and, according to Slonczewski-Berger 
formula [30, hig h switching current density. PMA 
in CoFeB/MgO [3] and in Co2FeAl/MgO MTJs .16] 
demonstrated recently with ultra-thin layers of CoFeB 
and Co2FeAl originates from the electrode/spacer inter¬ 
face, not the volume of the electrodes. Finally, TMR in 
MnaGa/MgO MTJs was found very small, far below the 
application range. 

In present letter we study BZ filtering mechanism 
of (exponentially) large TMR in MTJ based on hep 
Co(0001) electrodes and h-BN spacer. h-BN is an ul- 
trahigh chemically stable semiconductor with band gap 
of ~6 eV that has prefect lattice matching with hep Co. 
Growth of single and multiple layers of h-BN on Co(0001) 
has been recently demonstrated [0. To the best of 
our knowledge only few studies of Co/BN/Co s yst em ex¬ 
ist, mainly with single sheet of BN spacer |23l425| . The 
FM/BN/Gr/FM(111) MTJ has been suggested in |2ji], 
where the ferromagnet (FM) like fee Ni or Co do not 
have PMA . High TMR in this design is achieved by us¬ 
ing graphite (Gr) as a spin filter, while overall resistance 
of MTJ is regulated by thickness of h-BN. In present work 
we show that complex BN/Gr structure is not needed and 
h-BN itself can act as a spin filter to produce high TMR. 
Also, h-BN is more oxidation and intercalation resistant 
than graphite. 

We calculated transmission and I-V curves of the 
Co(0001)/h-BN/Co(0001) MTJ using nonequilibrium 
Green’s function (NEGF) approach [27] developed within 


the LDA-based LMTO-ASA formalism [28j,l29j. Geome¬ 
try of the system consist of Nbn h-BN sheets sandwiched 
between two semi-infinite hep Co electrodes. The in-plane 
lattice constant, 2.50 A, and distances between Co lay¬ 
ers, 2.03 A, and between h-BN layers, 3.33 A, were set to 
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FIG. 2: (color online). Left panels: transmission functions 
(a) T n , (b) Tqq, and (c) 7qq calculated at zero voltage for 
Co(0001)/h-BN/Co(0001) MTJ with N B n = 4, 6, 8, and 10. 
Right panels: currents 7qq, Iqq, and hi f° r (d) Nbn = 6 and 
(e) Nbn = 8, and (f) corresponding TMR shown as function 
of bias voltage. Inserts in (d) and (e) show I-V curves at 
smaller scale. Note logarithmic scale in (a), (b), (c) and (f). 


experimental values for bulk hep Co and h-BN. Positions 
of surface B and N atoms relative to surface Co layer 
were calculated by using VASP molecular dynamic pro¬ 
gram [30]. In lowest energy configuration N and B atoms 
occupy top and lrcp sites with distance between Co and 
BN layers equals to 3.31 A (see, e.g., Q for definition of 
top, fee and hep sites on the Co(0001) surface). Two con¬ 
figurations (N,B)=(top,hcp) and (N,B)=(top,fcc) have 
similar energies, in agreement with previous calculations 
[3 . 12H ] and experimental observations [22] . 

Transmission functions Tqq, Tqq, and Tqq calculated at 
zero bias voltage for Co(0001)/lr-BN/Co(0001) MTJ with 
Nbn = 4, 6, 8, and 10 are shown on left panels of Fig 
2. Here Tqq and Tqq are transmission in majority and 
minority channels calculated for parallel magnetization 
configuration of electrodes and Tqq is transmission cal¬ 
culated for antiparallcl configuration. The I-V curves of 
/qq, In, an d III currents per unit cell area are shown on 
Fig 2(d) and Fig 2(e) for Nbn = 6 and 8. Corresponding 
TMR = (Jqq + In ~ 2/qq)/2/qq is shown on Fig 2(f) as 
function of bias voltage, V. As seen on Fig 2(d),(e) the 
/qq at small voltages is about 7 times larger compared to 
/qq (and even more so for Jqq) resulting in ~ 250% TMR 
at V < 0.1 V [see Fig 2(f)]. This is mostly due to the fact 
that density of states (DOS) of minority electrons of hep 
Co at Fermi energy, Ef , is significantly larger than DOS 
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of majority electrons [see Fig 3(a) and Fig 3(b)]. 
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FIG. 3: (color online), (a) majority and (b) minority electrons 
DOS of bulk hep Co at Ef averaged over the wave number 
along c-axis, k z , shown as function of in-plane fc|| in 2D Bril- 
louin zone (BZ). Note absence of majority electron states near 
K and M points, (c) bands of bulk h-BN calculated by full 
potential (FP) LMTO method, (d) selected complex bands of 
h-BN calculated by Green’s function approach within LMTO- 
ASA [28i. |29j shown at three fc|| = Y,M and K as function 
of ( k z — k° z ) 2 , where k z is the wave number corresponding to 
the lowest energy of conduction band at given fc|| [k z = 0 for 
fc|| = r and M, and k z = 0.5 for k\\ = K], Negative ( k z — k z ) 2 
represent evanescent modes with k z (E) = k z + i'y(E). 

More interesting, however, is the sharp increase of /j_j_ 
and TMR that occurs at V ~ 0.4 V for Nbn = 6 and at 
V ss 0.3 V for N B n = 8 (TMR is as high as 4000% at 
0.4 V for Nbn = 8 ). Three features of transmission func¬ 
tions T|-|-, T pi., and are essential for understanding 
high TMR of the system predicted by LDA at V > 0.3 
V. These features are (1) sharp drop of transmission func¬ 
tions at energy near the valence band maximum (Ey) of 
the h-BN slab, (2) alignment of the Fermi energy, Ep, 
of Co inside the h-BN band gap very close to Ey (LDA 
predicts A Epy = Ep — Ey m 0.2 eV), and (3) sharp drop 
of Tjq, occurs at higher energy compared to energy where 
the sharp drop of T n and T-jq occurs [as seen on Fig 2 
T 44 is several orders of magnitude (!) larger than Tjq. or 
T-pf- in energy window Ep — 0.6 eV < E < Ep — 0.2 eV]. 
The last feature is a consequence of the BZ filtering as 
will be described below in discussion of Fig 4. 

Let us first explain the sharp drop of transmission func¬ 
tions shown on Fig 2(a)-(c) at energy near Ey. On Fig 
3(c) and Fig 3(d) we show real bands and selected com¬ 
plex bands of h-BN. In bulk h-BN the mode with smallest 
attenuation constant 7 (E) = Im(k z ) [shown by dashed 
green line in Fig 3(d)] is the mode with fc|| = 0 for all 
energies, E, inside the band gap except energies in close 
vicinity of valence band maximum, Ey, or conduction 
band minimum, Eq■ This evanescent mode continuously 
connects two real h-BN bands shown along A — T line on 
Fig 3(c): highly dispersive conduction band with energy 
Ey + 4.8 eV at T point and valence band with energy 
Ey — 4.1 eV at T point. As seen on Fig 3(d) maxi¬ 
mum of this 7 (E) occurs at E ~ Ey. This explains 
why transmission functions shown in Fig 2 have a ten¬ 
dency to smoothly decrease (on logarithmic scale) when 
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energy decreases from Ec to Ey with a minimum near 
Ey. On the other hand, since the valence band maximum 
(VBM) of h-BN corresponds to a state with momentum 
near the K point, for E slightly above Ey the evanes¬ 
cent mode with smallest attenuation constant [shown by 
dashed red line in Fig 3(d)] is the mode with fc|| close to 
the K point. Since the valence band along the H — K 
line with energy near Ey is almost flat [see Fig 3(c),(d)], 
corresponding 7 (E) increases very fast when E increases 
above Ey. Thus the behavior of the smallest attenua¬ 
tion constant of BN inside the band gap explains general 
shape of transmission functions of Co/h-BN/Co MTJ for 
E inside the h-BN band gap and their sharp drop near 
Ey. In other words, specific features of the complex band 
structure of BN described above are responsible for ultra¬ 
sensitivity of spin-dependent conductance and TMR of 
this MTJ to position of the Fermi energy inside the BN 
band gap. (Note that transmission functions of conven¬ 
tional MgO-based MTJs do not exhibit sharp features as 
function of energy since MgO has much simpler complex 
band structure as compared to that of h-BN.) 

Left panels of Fig 4 show bands of repeated slabs of 
9 hep Co and 5 h-BN layers calculated along symmetry 
lines of 2D BZ by LDA-based FP-LMTO method. Green 
and blue colors of bands on Fig 4 are mixed according 
to the values of projections of the wave functions to Co 
and BN atomic orbitals, correspondingly. Fig 4(a),(d) 
show that highest energy at which majority Co states 
and valence h-BN states co-exist at the same fcn is E' v = 
Ep — 0.6 eV, while highest energy at which minority Co 
states and valence h-BN states co-exist is the VBM of 
h-BN Ey = Ep — 0.2 eV. This explains why the sharp 
drop of T jq and T-jq (together with T^) shown on Fig 2 
occur at two different energies: Ep — 0.2 eV and Ep — 0.6 
eV, correspondingly. 

Let us stress again that it is the BZ filtering mecha¬ 
nism which is responsible for several orders of magnitude 
difference of T ^ transmission compared to 7)-j~ (and T-j-j,) 
at energies in the range (Ep — 0.2 eV, Ep — 0.6 eV). 
The real states of BN (or ’hot spot’) near the K point at 
these energies provide metallic transmission for minority 
Co electrons, X 44 , while absence of majority Co states in 
vicinity of K point at these energies makes T )-|- and Tjq 
exponentially suppressed. 

Our NEGF calculations show that for small bias volt¬ 
age, V, applied to Co/BN/Co system the equilibrium 
transmissions functions shown on Fig 2 do not change 
much if Ep is replaced by average (Ep) = (Ep + Ep)/ 2, 
where Ep and Ep are Fermi energies of the left and right 
electrodes ( Ep—Ep = eV). Thus, the onset of sharp rise 
of could be estimated at V « 2A Epy/e « 0.4 V [that 
agrees with actual NEGF onset at 0.3 — 0.4 V shown on 
Fig 2(d),(e)] and onset of sharp rise of Jp|- and Iy |- could 
be estimated at V ~ 2 (Ep — E' v )/e « 1.2 V, which is 
beyond the scale of voltages shown on Fig 2. 

Our results are very sensitive to the position of the 
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FIG. 4: (color online). Top panels: majority bands of repeated 
slabs of (a) 9Co/5BN and (b) 5Co/5BN calculated using the 
LDA theory, and (c) 5Co/5BN calculated using the QSGW 
theory (y-axis on all graphs is E — Ef in eV). Bottom panels: 
corresponding minority bands. 


Fermi energy relative to the bottom edge of the BN band 
gap, AEpv, since onset of high TMR occurs at volt¬ 
age V — 2AEpy/e, yet V cannot exceed the electri¬ 
cal breakdown limit. We studied dependence of A Epy 
on variations in geometry of the systems by performing 
FP-LMTO calculations within LDA for repeated slabs of 
Nbn h-BN and Nco hep Co layers (varying Nbn and 
Nco from 4 to 12) for all configurations of B and N 
atoms (B,N)=(top,hcp), (top,fee), (hep,top), (hep,fee), 
(fee,top), and (fee,hep) and for reasonable variations of 
positions of surface B, N and Co atoms. The value 
of A Epv was found to be stable within LDA/DFT, 
A Epy « 0.2 eV, in all above cases for sufficiently thick 
Co slab with Nc 0 > 7. 

On the other hand, it is well known that LDA underes¬ 
timates band gaps of semiconductors and poorly describes 
the Schottky barrier heights 31]. To evaluate A Epy on 
the level beyond LDA we performed calculations for pe¬ 
riodic slab of 5 hep Co and 5 h-BN layers by using the 
QSGW theory that is known to describes band gaps and 
other properties of materials with moderate e-e correla¬ 
tions significantly better than LDA fjil - ljjt . (The restric¬ 
tion to 5Co/5BN slab size was due to heavy computa¬ 
tional costs of QSGW.) Fig 4 shows the band structure 
of majority and minority electrons calculated for periodic 
9Co/5BN and 5Co/5BN slabs within the LDA theory, 
and for 5Co/5BN slabs within the QSGW theory. It is 
seen that reducing the size of Co slab from 9 to 5 lay¬ 
ers leads to moderate increase of LDA A Epy from 0.2 
eV to 0.5 eV (as mentioned above A Epy converges to 
ss 0.2 eV starting with Nco > 7), whereas more accu¬ 
rate inclusion of the e-e correlations within the QSGW 
theory results in dramatic increase of A Epy to as large 
as 2.2 eV. Thus, taking into account correlations on be¬ 
yond LDA level is very important for accurate prediction 
of the Fermi energy alignment (A Epy) at Co/li-BN in¬ 


terface and, consequently, for designing the Co/h-BN/Co 
MTJ with high TMR. 

In order to decrease A Epy from 2.2 eV predicted by 
QSGW to practical levels one can p-dope the h-BN (e.g. 
by Mg). It was recently shown that Mg acceptor level 
in Mg-doped h-BN (h-BN:Mg) is as low as 0.031 eV [35}, 
so varying the Mg concentration could bring AEpy to 
desired range of AEpy < 0.1 eV. Note that for ideal 
alignment of AEpy = 0 the TMR could be as high as 
several orders of magnitude at very low voltages [see Fig 
2 ]- 


In conclusion, general ’Brillouin zone spin filtering’ 
mechanism of enhanced TMR is described for magnetic 
tunnel junctions and studied on an example of MTJ with 
hep Co electrodes and h-BN spacer. We suggest new MTJ 
for STT-MRAM and STO applications based on hep Co- 
Pt or Co-Pd disordered alloy electrodes and p-doped h- 
BN spacer. PMA of such low symmetry hep electrodes 
originates from the whole volume rather than interfaces 
and does not require chemical ordering at atomic scale 
promoting small device-to-device variations, a property 
important for high yield and low cost. Concentration 
of heavy Pt or Pd atoms could be balanced between re¬ 
quirements of strong PMA and small damping constant 
of the alloy (damping constant of pure hep Co is small, 
a hcpCo = 0.004 [21 at 300K). 

Owning to the specific complex band structure of the 
h-BN the spin-dependent tunneling conductance of the 
system is ultra-sensitive to small variations of the Fermi 
energy position inside the BN band gap. Our LDA-based 
NEGF calculations for Co(0001)/h-BN/Co(0001) MTJ 
shows high TMR at V > 0.3 V due to Brillouin zone 
filtering mechanism. Critical property needed for high 
TMR in this MTJ is an alignment of Co Ep very close 
to VBM of h-BN. LDA predicts AEpy ss 0.2 eV, while 
AEpy calculated by more accurate QSGW theory in¬ 
creases to 2.2 eV. Thus p-doping (e.g. by Mg) of the 
h-BN is needed to reduce AEpy to practical range of be¬ 
low 0.1 eV to ensure high TMR (that can be as high as 
several orders of magnitude) at low voltages. By varying 
the Mg concentration (and thus A Epy) and h-BN slab 
thickness one can independently regulate the TMR and 
overall resistance of the MTJ. 
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